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We consider the mixing of muon and tau neutrinos with sterile fermion fields 
, propagating in extra dimensions in the context of core collapse supernova 

\^ \ physics, extending the analysis of the electron neutrino case done in a previous 

■ work. We show that the potentially dramatic modifications to the supernova 
. evolution are prevented by a mechanism of feedback, so that no severe bounds 

^) \ on the parameters of the extra dimensions need to be imposed. Nevertheless, 

■ the supernova core evolution is significantly modified. We discuss the conse- 
quences on the delayed explosion mechanism and the compatibility with the 



SN 1987A signal. Then, for the cases of both u^^^ and Ve mixing with bulk 
fermions, we analyse the distinctive features of the signal on Earth. 



]^ ■ 1 Introduction 

H . . m 

I The study of sterile neutrinos propagating in the bulk of extra space dimensions [T] and 

of core collapse supernova evolution [2] establishes a nice link that can have interesting 
implications on both sides, even though the role of sterile neutrinos is not believed to be 
primary in solar and atmospheric contexts 1^. We consider indeed a framework [3] 
in which the solar and atmospheric oscillation signals are independently accounted for 
in the standard way. In addition, a small mixing with the 4-dimensional Kaluza-Klein 
(KK) states (the 4d correspondent of the extra dimensional field) is the source of the 
non-standard effects on supernova (SN) physics, since it implies resonant oscillations with 
ordinary Standard Model neutrinos via matter effects. The point is that we have a large 
number of resonances that can enhance the effect, even if the mixing angles are small. As 
KK modes do not interact with ordinary matter and are therefore not trapped, conversions 
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lead to a faster cooling and deleptonization of the SN dense core. The observed SN 1987A 
signal [S] requires such an energy loss not to take place with a time scale much shorter than 
10 sec. However, the consequent potentially severe bounds on the extra dimension size [7] 
are widely relaxed by a feedback mechanism jSl IH] that self-limits the conversion. In a 
previous work jH], we analysed the possibility of a mixing between the electron neutrinos 
and a sterile bulk fermion and we explicitly showed that the mechanism indeed prevents 
the potentially dramatic cooling in a wide portion of the parameter space. We also pointed 
out that new interesting effects on the SN dynamics can arise in the enlarged parameter 
space. We speculated that the enrichment of the Vf. component in the outgoing neutrino 
flux might help the delayed shock mechanism and drive the explosion.^ Furthermore, we 
found that while this picture is in good agreement with the observed SN 1987A signal, the 
flux composition can be significantly modified. 

In the first part of the present work we address the case of muon and tau neutrinos. 
In the standard picture of collapse and cooling of the SN core, their role is rather different 
from the electron ones. In fact, while electron neutrinos are copiously produced by the 
neutronization process during collapse and build up a degenerate Vf. sea with large chemical 
potential /i^^ (therefore electron antineutrinos are a negligible fraction in these conditions), 
muon and tau neutrinos and antineutrinos are produced in pairs and have almost vanishing 
chemical potential. Moreover, lacking in charged current interactions, they have a faster 
diffusion and give an important contribution to the transport of thermal energy out of the 
core. Also, their lepton number does not play a significant role in the diffusion dynamics. 

The mixing of tau or muon neutrinos with bulk fermions opens an unconventional 
escape channel, which can affect the SN core cooling and deleptonization in a peculiar 
way. We will see that the usual energy loss bounds can be relaxed, as in the electron 
case, by the action of a feedback. We will then be interested in the phenomenological 
consequences on the revival of the shock and in the implications for the SN 1987A signal. 

In the last part of the paper we discuss the peculiar signatures of the considered effects 
on the neutrino signal on the Earth. This we do not only for the case of f^, but also 
for the case of v^. conversion considered in [2]. 

In the following we will assume that the sterile fermions mix predominantly with the 
tau neutrino. This is the case if the SM/bulk neutrino couplings reflect the hierarchical 
structure of the SM Yukawa couplings. Actually, the inclusion of a significant muon neu- 
trino mixing with bulk fields would just unnecessarily complicate the core dynamics. The 
point is that the oscillations into bulk neutrinos build up a non-zero tau (or muon) lepton 
number. The tau neutrino abundance is simply determined by the tau lepton number, since 
tau leptons are not thermally produced, as their mass is large with respect to the energies 
involved (T ~ 30MeV). On the contrary, the muon neutrino abundance would also de- 
pend on the non-negligible muon fraction through the /3-equilibrium processes p <-> n 
and n •f-^ pz/^. One can also wonder whether the Vpr'^^T mixing could reintroduce the 
muon fraction in the problem. This is not the case, since that mixing is largely suppressed 
by the difference at one loop of the muon and tau matter potential ^0] and by an even 
larger difference originating from the new effect. In fact, as soon as the resonant conver- 
sion into the bulk begins, a non-vanishing tau lepton number is generated and thus the 

^Thc difficulty to obtain an explosion in simulations, however, can be explained by the complexity of 
the system to be simulated, so that no new physics is in principle necessary. 
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Figure 1: Excluded regions in the m-R plane, for the case of Vr (i^ii) escape into the 
hulk. The red dashed line (corresponding to m^R = 10~^^ e\ ) marks the lower border of 
the region where the supernova core evolution is sensibly modified. In the region on the 
right (1/R> IkeVj the effect is not active: neutrinos do not cross a sufficient number of 
resonances. 

matter interactions produce different effective masses for tau and muon (anti)neutrinos. 
For sufficiently large mixing with the extra dimension fields, as is our case, the oscillations 
are suppressed before they significantly influence the neutrino fractions. Notwithstanding 
the complications described above, we do not expect a qualitatively different behaviour in 
the presence of muons; in section 0] we will thus deduce the case of mixing with bulk 
fermions from our results. 

The paper is organized as follows. In section [21 we briefly describe the general extra 
dimensional scenario and set the parameter space in which we are interested. In section |21 
we first sketch the expected effect of such new physics on the evolution and then present 
the numerical model that we use. Finally, in section |3] we address the new phenomenology 
based on the results of our calculation and we compare the expected neutrino flux on 
Earth in the cases of fg, and I'r conversions. In section El we sum up. 

2 Framework 

The cornerstones of our large extra dimension scenario are the following. We consider a 
single gravitational extra dimension ^, probed by a sterile fermionic field, which results 
from a 4d point of view as a KK tower of states, predominantly mixed with tau neutrinos. 

The parameter space we are interested in is illustrated in fig. ^ The radius R of the 
extra dimension, which sets the scale of the KK spectrum, is allowed to lie in the broad 

^More extra dimensions could be present at higher energies. This scenario is compatible both with a 
fundamental gravity scale at ---^ 1 TeV |11| and with astrophysical and cosmological bounds |12|. 
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range 10"'^ eV ^ ^ 1 keV. The lower bound comes from direct gravity probes, the upper 
one from requiring that a sufficient number of resonances be crossed by neutrinos travehing 
in the core, a condition for observable effects to arise. The mixing angle of the Uj- with 
the A;-th eigenstate of mass is generically parametrized as 6k — /"^ , which defines 
the parameter m; we only assume that the density of KK states is proportional to R, a 
feature that is quite model- independent. The condition mR <C 1 ensures the smallness of 
the mixings, while keeping subleading effects under control [HI requires mR <C 2 x 10~^. 
The above requests identify the unshaded region in fig. ^ In this region the new effects 
enter the analysis only through the crucial parameter rn^R. 

If feedback effects were not taken into account, only values of m^R up to 10~^^eV 
would be compatible with the energy loss argument. For example, for m^R = 5 x 10"^ eV, 
the core would cool down in about 10~^ sec. On the contrary, including the feedback, any 
value of m^R up to 5 x 10~^ eV turns out to be acceptable, while beyond that point the 
core is cooled by otherwise subleading effects. Thus, a new portion of the parameter space, 
10~^'^eV < nn?R ^ 5 x 10~^ eV, is opened once the feedback is taken into account. ^ In 
the following we will focus on the range 

10"^^ eV < m^ii < 5 X 10-^ eV , (1) 

where the neutrino escape affects the core evolution in a sizeable way. 

Tau neutrinos in the SN core experience a Mikheyev-Smirnov-Wolfenstein (MSW) 
potential given by 

Vr = ^GyUb Q^e + 2y., + - ^) , (2) 

where Gf is the Fermi constant, is the baryon number density, and Yx is the net 
number fraction per baryon of the species x: Yx = {Nx — Nx)/nB- It will also be useful 
to consider the lepton number fractions Y^^^^^^ = Ye^f_i^r + Y^^^^,^^^^. 

At the beginning Yg is the dominant term and, since it is not expected to be larger 
than 0.4 in the inner core [2j, is clearly negative. The tau antineutrinos acquire an 
effective squared mass m^jj = 2EpV, where Ep is the antineutrino energy. The effective 
mass changes slightly along the neutrino path: whenever m^s equals the mass of one of the 
KK sterile states, a resonance occurs and the tau antineutrino has a certain probability 
of oscillating into a bulk state and escaping from the core. As a result, the tau lepton 
number increases and the potential becomes less negative. If the potential were positive 
the resonances would be met by neutrinos and the potential would get smaller. Note also 
that Y^^ is zero at the beginning but a positive value will be generated by the Ur escape 
itself. This term, which has a different weight for the f^, is responsible for the suppression 
of the vacuum v^-Vr oscillations. 

The disappearance probability of a Vr into the bulk when travelling a distance L is 
given by the usual Landau-Zener formula 

P(.,^bum=.L^„'ij(My'\ (3) 



^Compared with the electron neutrino escape case, the overall shift to higher values is justified (i) by the 
smaller number density of Vt with respect to u^, since they are thermal, not degenerate, at the beginning; 
(ii) marginally, by the fact that Vr diffusion is faster. 
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As is apparent, the parameter Tm?R sets the overah magnitude of the escape effect. 

3 Model and evolution 

Having introduced all the ingredients of the problem, let us first discuss the basic features of 
the expected evolution, compared when necessary with the case of Ue conversion addressed 
in 0. 

As already mentioned, the matter MSW potential is initially negative in the whole 
inner core so that only Uj- cross resonances. As conversions push Vr towards zero and 
deleptonization pulls it down to its negative minimum value, a non-trivial interplay be- 
tween diffusion and the new effect takes place. Nevertheless, the main features of the 
evolution can be easily understood. First of all, since the new physics is in the u^- sector, 
the Yl^ evolution is only indirectly affected, mainly via the modifications to the temper- 
ature. As the VeS are almost degenerate, we expect their diffusion to be only slightly 
changed and in particular the deleptonization time scale to be unaffected. On the other 
hand, the U-j- escape generates a positive Y^,^, the balance Vt-Vt is broken and a positive 
chemical potential ^i^^ arises. This in turn inhibits the escape itself, both because Vr 
is lifted towards zero by the term Yy^ in eq. ((2j and, more important, because the v^- 
abundance is suppressed in the presence of the chemical potential. In particular, has 
no chances to switch to positive values and only do convert during all the evolution. 
Now, how do we expect the T evolution to be affected? In a first fast phase, as long as 
the conversion is dominant, a portion of the thermal energy is absorbed both by bulk neu- 
trinos and by the sea of increasingly degenerate z^t-'s. During all the following evolution, 
the thermal energy that the z^t-'s escape keeps on transferring to tau neutrinos is carried 
outside by their diffusion. Thus, we expect the T evolution to be sped up with respect to 
the standard case. Moreover, we expect the escape cooling channel to be more effective 
than in [H] as the Vt diffusion time scale is shorter than for fg's. Finally, the tau lepton 
number that arises in the inner core diffuses out into the mantle and may significantly 
influence the outgoing flux composition. 

To study the above picture in some more detail, we use a simplified model that includes 
all the relevant physics, although it is not intended as a complete description of the SN 
core dynamics. A more detailed analysis would require a complete simulation of the core 
and of the external layers, with the inclusion of the new effect. 

We follow all the core evolution in terms of the temperature T, the electron leptonic 
fraction 1^^ and the tau leptonic fraction Yi^^ (that is equal to the neutrino fraction 
Yy^\ One can express all the other relevant observables (chemical potentials, entropy, 
fluxes, . . . ) simply as functions of these quantities, using thermodynamic equilibrium, beta 
equilibrium and charge neutrality. As in [Hj, we consider a typical SN core characterized 
by a mass of 1.5 Mq and a radius of 12.7 km. For the matter density, we use a static 
profile p{r\ justified by the relative hydrodynamical stability of the inner core in which 
we are interested. Explicitly p{r^ = Pc/{^ + {i'/^)'^)-, with pc = 7.5 x 10^^ g/cm^ and 
4/3 TT T Pq — 1.1 • As for the equation of state, we consider a core made only of ?i, p, 
e^, i^e,fj.,T, ^e,ii,T and 7. The effective nucleon mass is expressed as = jj^rp^^^jj-^, where 
rriN = 939 MeV is the vacuum value, /3o is chosen to be 0.5 and po = 3 x 10^^ g/cva? is the 
reference nuclear density. As for the mean free paths X^^ and Xy^^^y^, in general they would 
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exhibit a non-trivial dependence on the neutrino energy E^, and the evolution variables (p, 
T, Yl^, besides on several other aspects related to the medium in which the neutrinos 
diffuse. However, the essential features can be grasped by assuming a simple inverse 
quadratic dependence on the neutrino energy for all species and incorporating an inverse 
dependence on matter density for muon and tau neutrinos, while keeping constant with 
density the mean free path of electron neutrinos. We checked these simplified assumptions 
versus the more complete modellings of JHl- Moreover, the above choices allow us to 
obtain an evolution whose main features and timescales agree with the results of more 
sophisticated analyses. The expressions we use are the following:^ 

E"^ P E^ 

XueiEu) = , ^v^,vAEy,r) = , (4) 

taking = 1.2 cm and = 2.8 cm at the reference energy E^^q = 260 MeV and 

reference density pc = 7.5 x 10^^ g/civ? . The initial profiles for T and Yi^ are typical ones 
as taken in Of course, the profile for Y^^ is zero in all the core at the beginning. 

Finally, with all the above prescriptions (see [0] for a more detailed discussion), the 
evolution equations for Y^^, Y^^ and T read 

OYl 

^B—Q^ = V • (oe Vpue) (5a) 
= V • (a. V^.J + -^m''R^\ T^/^Fy, f-^) (5b) 



ubT— = ae [Vpu.j +ar[yiJ-urJ + V • ( (oe + + a^)— VT' 



(5c) 



where Oj = A5^^£'^Q/(67r^), i = e,p,T; s is the entropy per baryon and the Fermi integrals 
are defined as Fq{y) = dx-^^^^ . 

One can recognize the standard transport equation for electron neutrinos in eq. (|SH)i 
while eq. (jSh)) for the tau neutrinos shows the additional escape term that comes from the 
integration over neutrino energies of the probability in eq. Q, for L equal to a mean free 

path. The first line of the temperature eq. ((Ht) features the heating terms Oj Puj^ 
{j = e, t) associated to the degradation of the degeneracy energy of neutrinos reaching 
regions with lower chemical potential. The novelty is the presence of such a term for 
tau neutrinos also, since they possess a non-vanishing chemical potential now. After the 
standard temperature diff'usion term V • (VT^), last comes the contribution from tau 
antineutrinos escape, which is a cooling term: its first part is related to the energy fiux 
into bulk when an antineutrino is lost, the second one is associated with the (negative) tau 
lepton number fiux into the bulk and thus with the storing of energy in the tau neutrino 
degenerate sea. 

The numerical solutions of the above equations, for several choices of the parameter 
m?R in the range of eq. confirm the expectations presented at the beginning of this 

*With respect to the electron neutrino case of [3| we refined the value of A[J^ and we introduced the 
density dependence for muon and tau neutrinos. 
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Figure 2: Profiles of the MSW potential for tau neutrinos. The thick line is the initial 
configuration corresponding to the typical initial profiles assumed in f^; the thin lines are 
snapshots at 1, 5, 10 and 20 sees, following the arrows, for the case m?R = 10"'^ eV . 

section. Essentially, the electron neutrino density follows the standard evolution while the 
temperature fall is faster. The tau leptonic fraction Y^^ rapidly grows from zero up to a 
certain profile (the higher the larger m?R, but always below about 0.1) and then lowers 
on typical diffusion timescales. In parallel, the chemical potential fiy^ quickly grows up to 
values as large as those of fiy^ (of order 200 MeV) and then clears out with diffusion. A 
typical evolution of the matter potential is shown in fig. |2l at the beginning it is pushed 
towards zero by the rising of Y^^ but, the initial potential being well below zero, the 
available energy is not sufficient for a complete zeroing, and diffusion soon reverses the 
trend. 

4 Phenomenology 

Figures El and show the energy and the tau lepton number drained out of the inner core 
in the first ten seconds of evolution, separating the portion that is actually radiated in the 
visible neutrino flux and the one that is irreparably lost into the invisible channel provided 
by the bulk. The first ten seconds are the most interesting interval since this is the lapse 
of time during which the relevant neutrino signal is produced. 

As for the energy, one sees that a sizeable portion can be lost into the invisible chan- 
nel, especially for high values of rn^R. At the same time, the total amount of drained 
energy increases significantly with nn?R, as a consequence of the faster cooling of the core. 
However, balancing the two effects, what is important is that the reduction of the portion 
emitted into the visible channel is always limited to ~ 20% of the standard case, which is 
well acceptable in the light of the need of accounting for the SN 1987A signal: no direct 
upper bound on the parameter m^R then needs to be imposed. 

The other significant result is the emission of a net tau lepton number, counterbalancing 
the tau antineutrinos that escape into the bulk. The emission can grow up to values ~ 10 
times larger than the corresponding fiux, for m?R at the end of the considered range. 
As for its time distribution, fig. El shows that the v^- emission is concentrated in the very 



7 



^Q-13^Q-12^Q-ll^Q-10^Q-9 ^^-6 

m^R (eV) 



^Q-13^Q-12^Q-ll^Q-10^Q-9 ^^-6 

m^^ (eV) 



Figure 3: Energy drained out of the SN 
core in the first 10 sees. The total energy 
(red solid line) is shown together with its 
visible (upper blue dashed line) and invisi- 
ble (lower blue dashed line) components as 
a function of m^R. 



Figure 4: Tau lepton number emission 
in the visible channel (blue dashed line) in 
the first 10 sees. The electron lepton num- 
ber emission (green dot-dashed line) is also 
shown for comparison. 
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Figure 5: The tau lepton number flux versus time during the first 5 sees, for increasing 
values of m^R = 10"^^^' ^' ^\ from bottom to top. The normalization is given by the 
value of the electron lepton number flux in absence of new physics. 

first seconds. An analogous behaviour is exhibited by the total energy emission, so that we 
can foresee a peculiar time dependence of the final neutrino number flux on Earth, more 
peaked at earlier times. Letting aside this general expectation, we focus in the following 
on time- integrated quantities. 

Taking our outputs of the modified core evolution as the starting point, we now want to 
outline the subsequent vicissitudes and the ultimate fate of the neutrinos from the cooling 
phase in order to derive the possible consequences both for supernova physics and for the 
detectable signal on Earth. 

First is the issue of how the energy flux and the sizeable I'r net flux emitted from 
the inner core are redistributed by the subsequent diffusion that takes place up to the 



neutrino sphere (after which the neutrinos stream freely) . A detailed examination of this 
phase would require studying the evolution of the mantle also; see e.g. Jl]. What we 
do is to assume that the energy E emitted from the inner core ends almost equiparted 
in neutrinos and antineutrinos of all flavours, once the portion carried by the lepton 
number excesses has been subtracted. On the other hand, the (positive) lepton num- 
bers Lf> and are separately conserved in the reprocessing. With these ingredients 
one can then determine the neutrino and antineutrino number fluxes in each flavour: 
np, ^ ^(^-Ei(^!^,)^i)/(3 + Ei(^!^,)/(^pJ)> - Lj + np^, where i,j = e,fi,T. 
For the average energies of the different families we use the typical values (E^^) « 13MeV, 
{EpJ Ki 16MeV and « 23MeV.^ It has recently been pointed out [K] that the 

mean energies of muon and tau neutrino and antineutrino can be closer to the electron 
antineutrino's, at the level of less than 1.2 {Ep^) during the core cooling phase. Neverthe- 
less, such a variation does not change our analysis. We will comment on this issue later 
on. 

Thus, coming out from the neutrino spheres, the flux displays a peculiar composition 
(see fig. El): for large values of m^R, the i^^'s are ~ 5 times more than in the case without 
new physics, at the expense of the i^e's and of all flavours of antineutrinos (all reduced by 
up to ~ 60%). These reductions, although not so large, go into the direction opposite to 
the need of rejuvenating the stalling shock wave and thus helping the explosion, since muon 
and tau neutrinos are much less efficient than the electron ones in transferring their energy 
to the lingering matter. On the other hand, this qualitative argument does not imply a 
strong drawback. The actual explosion mechanism may be hidden in the complexity of 
the system, whose simulation is still a non-trivial task. 

Next, neutrinos and antineutrinos undergo the matter flavour oscillation in the periph- 
eric low-density region of the star and the vacuum flavour oscillations in the journey from 
the supernova to Earth. We illustrate the outcome in the case of normal neutrino mass 
hierarchy. As for the oscillation parameters, we use 623 — 45° [1] and 9i2 — 32° while 
we consider values for ^13 compatible with the Chooz bound J7]- The first transition at 
very high densities is produced by the small matter potential difference between tau and 
muon neutrinos, which is due to the mass difference of the corresponding leptons. This 

transition is found to be completely adiabatic and switches the ^t'*^, 'i^V flavour eigenstates 

to the states rotated by the angle 623, labelled ''u and At lower densities, two more 
resonances are met, due to the charged current contribution to the electron matter po- 
tential. In the case of normal hierarchy, these resonances are both active in the neutrino 
sector. The Ve-v^i resonance turns out to be adiabatic under our assumptions while for the 
Vg-Vr' we consider a level crossing probability Ph, which can vary from to 1 depending 
on the value of ^13. For a more detailed analysis, we refer the reader to JH]. Finally, 
neglecting terms proportional to sin ^13, the neutrino and antineutrino fluxes F reaching 

^We are assuming that the average energies of the neutrinos at the neutrino sphere are not affected by 
the new physics effect, see 
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Figure 6: Indicative composition of the neutrino flux leaving the neutrino spheres as a 
function of m?R. 



the Earth's surface are related to the neutrino sphere fluxes F by the simple relations 

= PhsI^F^, + cI^F;^ + (1 - PH)s'i2F^ (6a) 
= {PhcIAs + (1 - Ph)sIs)F^ + sl,cl,Fli + {PhsI, + (1 - Ph)cI,cI,)F^ (6b) 
Fr = {PHci2sl3 + (1 - Ph)42)F^ + s?24^° + (PhcI^ + (1 - PH)ci,sl,)F^ (6c) 



Fe = c22F° + 4FP (7a) 

Fp = ^?24^° + ^L^^" + c?24^° (7b) 
= 42sl3F^ + clsFll + cI,sI,f9 (7c) 

where Sij and Cjj are the sine and cosine of the angle 9ij . 

The indicative final composition of the flux reaching the surface of the Earth is il- 
lustrated in figs. and [7f for the boundary values Ph = 1 and Ph = 0, corresponding 
respectively to sin^ ^13 < 10"^ and sin^ ^13 > 10~^. In the first case, the strong u-j- ex- 
cess is distributed between muon and tau neutrinos, while the electron fraction decreases. 
The situation is different in the adiabatic {Ph = 0) case, where we expect an higher Ve 
component. A general important feature is the z^e flux reduction that can amount to up 
to ~ 60%. Given the limited statistics of the 1987 event, the overall uncertainties on the 
expected neutrino fluxes (as predicted by full simulations, also based on assumptions on 
the progenitor star) and finally the approximate nature of the supernova core evolution 
adopted here, such a reduction is still compatible with present observations. Nevertheless, 
this feature is an interesting and potentially challenging one, and deserves to be addressed 
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Figure 7: Indicative composition of the neutrino flux reaching the Earth's surface, as 
a function of rn^R, for the different cases of escape into the bulk: Ve (first column), 
(second column), (third column); for Ph = 1 (first line) and for Ph = (second line). 

more closely in case of a future SN event with higher statistics. A similar remark also 
holds for the time structure of the neutrino signal mentioned above. 

Our analysis has been limited so far to the case of tau neutrinos, but can be readily 
extended to the case of mixing of muon neutrinos with the bulk fermions. Indeed, no 
relevant modifications are needed to the model of evolution described in sec. |31 so that the 
general conclusions reached for tau neutrinos are still valid. In particular, the indicative 
fluxes coming out of the neutrino spheres are given in first approximation by a plot similar 
to the one in fig. El with the excess in the channel. In turn, using eqs. © and eqs. Q, 
one can determine the compositions reaching the Earth's surface shown in figs. da andE^; 
the value of Ph is only marginally important in this case. 

Finally, we can recall the results of [HI for the case of the electron neutrinos oscillating 
into the bulk and derive the compositions in figs. [7^ and[7|l following the same steps. 

Based on the above results, one can easily point out the most important observable 
modifications of the expected flux in a detector. First of all, as the net effect is always 
an increase in the total emitted lepton number, a crucial signature is the enhanced ratio 
between neutrinos and antineutrinos. This is a general feature of supernova neutrino 
oscillations into extra dimensions that does not depend on which flavour has a significant 
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mixing with bulk neutrinos. What is case-dependent is which channel the enhancement 
shows up in. First one can focus on the electron channel, which will have high statistics 
in the future supernova events. Indeed, Ue are the dominant signal in the present and 
in most of the future detectors, while fe can be efficiently collected by large Cerenkov 
detectors, already existing (SK, SNO) or proposed (UNO), and liquid-argon detectors 
(Icarus, LANNDD) For several cases in fig.[7|the I'e/^e ratio is significantly enlarged. 
However, it is also clear from fig. [3 that, depending on the value of Ph (i.e. of ^13), the 
effect of z/g or Vr conversion into the bulk could turn out to be irrelevant for the fg/z^e 
channel in some cases, so that it would be necessary to measure the or v-r fluxes. This 
is a harder task, which is however feasible at SNO, future detectors (UNO, LANNDD, 
OMNIS) or even at large scintillator detectors (KamLAND, Borexino [SOI)- On the 
other hand, once an anomalous neutrino/antineutrino ratio is measured, how is it possible 
to distinguish which flavour does mix with the bulk? The answer lies in the energy 
spectrum. In fact, since the electron neutrinos are less energetic than the muon and 
tau ones when they leave the neutrino spheres, the effect of the Ve ^ ^^,t oscillation is to 
harden the spectrum on the Earth. The conversion into the bulk of the electron flavour, 
increasing the fg component flowing in the SN mantle, reduces this hardness, while the 
muon or tau conversion enhances it. 

We have illustrated our results for a specific point in the parameter space of the 
standard neutrino oscillations. The residual dependence on 9i2 (and ^23) is mild. On 
the other hand, the distribution of the neutrino flux enhancement in the three flavours 
depends significantly on the neutrino mass pattern. Disentangling the flavour structure of 
the mixing with bulk neutrinos would therefore require the knowledge of sign(Am23) j21j . 

As already mentioned, another source of uncertainty is the precise value of the mean 
energies of the neutrinos coming out from their neutrino spheres. We checked that the 
adoption of a {Eu^^T,y^,T) closer to {EpJ, as suggested in yields only minor modifi- 
cations of our result. Namely, each column in fig. IHl and figs, d is affected at most by 
a variation of 5 — 10% in composition, within the uncertainties of our simplified model. 
Thus, the analysis is unchanged and the general trend confirmed: we only remark that 
the closer the mean energies, the smaller the differences among the channels, although 
preserving the features described above. 

Finally, in order to predict the actual signal in the detectors, one should also consider 
matter-induced oscillations inside the Earth. This effect has the power of dramatically 
changing the expected signal, but of course depends on the geographical position of the 
detectors at the time of arrival of the signal. For a given path in the Earth, it is not 
difficult to include the matter effect (see |221 118j ) . 

While we have so far considered the case in which the neutrino escape is dominated 
by the oscillations in the bulk of one of the three neutrinos, it is also possible that all 
SM neutrinos are involved in the energy loss. In this case, the neutrino flux on the Earth 
would still be enhanced over the antineutrino one. However, the neutrino flux composition 
would also be determined by the relative size of the three mixing parameters mg, m^, rrir- 
The evolutions of the three neutrino abundances in the core would be coupled through 
the temperature and the matter potential. The feedback mechanisms would however still 
take place, but we expect the rrii-R parameter space for the three channels to be slightly 
reduced. 
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5 Conclusions 



We dealt with the mixing of the Standard Model neutrinos with bulk sterile fermion fields, 
a natural and general feature of models with extra dimensions. In particular, we focused 
on the muon and tau case, while the electron case was addressed in [^. This mixing 
provides an unconventional escape channel for neutrinos in the supernova core during 
the cooling phase, which could in principle give strong bounds on the parameters of the 
extra dimensions when facing the observed SN 1987A signal. In order to simplify the core 
dynamics, we restricted the analysis to the case of i^r, but we easily completed it to the 
case of i^^. 

Given the sign of the matter potential, the process consists in the escape of tau an- 
tineutrinos into the bulk, acting in parallel with diffusion. We showed that a feedback 
mechanism turns on, preventing extreme and unacceptable modifications of the evolu- 
tion to occur (in particular, the characteristic time scale of deleptonization and cooling, 
of order 10 sec, is preserved). The direct bound on the parameters of extra dimensions 
can then be relaxed by about 4 orders of magnitude. The upper limit on the relevant 
combination m^R is lifted from 10~^^eV to 5 x 10~^eV, the new limit now coming from 
subleading effects. In fig.Q we showed the allowed region in the m-R plane. The feedback 
is mainly due to the building up of a degenerate sea of trapped i^r and thus of a positive 
tau neutrino chemical potential that suppresses the number of v'r^s candidate to escape. 
For comparison, in the z^e (and Ue) case in jH] this kind of feedback acts together with the 
feedback on the matter potential in a sizeable portion of the core. 

Although in such a 'safe' way, the supernova core evolution is significantly affected 
by the escape process, and we studied the main modifications. First of all, the cooling is 
sped up with respect to the standard case. Secondly, a sizeable net tau lepton number is 
emitted by the core, concentrated in the very first seconds. Therefore a neutrino number 
flux more peaked at earlier times can be predicted. 

Next we speculated on the phenomenological consequences on the revival of the shock 
and on the neutrino signal reaching the Earth. As for the first point, since the flux 
is lowered by the reshuffling of the neutrino fluxes (although not largely), we suggested 
that the v^, Vr escape may reduce the energy transfer to matter and thus the push to the 
explosion. As for the second point, the compatibility with the SN 1987A signal is preserved 
once one is willing to accept a certain reduction of the component. Collecting the results 
for the fe, and escape cases, we indicated as signatures, in cases of future galactic 
SN events, a general dominance of neutrino flux over antineutrinos and, depending on the 
flavour that mixes with bulk flelds, some peculiar structures of the relative enhancements 
in the neutrino fluxes and a harder or softer z/g spectrum. 
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